arXiv: 1507.00453v2 [cond-mat.stat-mech] llJul2015 


Universal Critical Wrapping Probabilities 
in the Canonical Ensemble 


Hao Hu"-'’, Youjin Deng^ *’-* 

^National Laboratory for Physical Sciences at Microscale and Department of Modern Physics, 
University of Science and Technology of China, Hefei, Anhui 230026, China 
^State Key Laboratory of Theoretical Physics, Institute of Theoretical Physics, 

Chinese Academy of Sciences, Beijing 100190, China 


Abstract 

Universal dimensionless quantities, such as Binder ratios and wrapping probabilities, play an 
important role in the study of critical phenomena. We study the finite-size scaling behavior 
of the wrapping probability for the Potts model in the random-cluster representation, under the 
constraint that the total number of occupied bonds is fixed, so that the canonical ensemble applies. 
We derive that, in the limit L — > oo, the critical values of the wrapping probability are different 
from those of the unconstrained model, i.e. the model in the grand-canonical ensemble, but still 
universal, for systems with 2y, - d > Q where y, - 1/v is the thermal renormalization exponent 
and d is the spatial dimension. Similar modifications apply to other dimensionless quantities, 
such as Binder ratios. For systems with 2y, -d <0, these quantities share same critical universal 
values in the two ensembles. It is also derived that new finite-size corrections are induced. These 
findings apply more generally to systems in the canonical ensemble, e.g. the dilute Potts model 
with a fixed total number of vacancies. Finally, we formulate an efficient cluster-type algorithm 
for the canonical ensemble, and confirm these predictions by extensive simulations. 

Keywords: Universal dimensionless quantity; Canonical ensemble; Finite-size scaling; Potts 
model 
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1. Introduction 

The fact that critical systems tend to display asymptotic scale invariance relates to the exis¬ 
tence of nontrivial dimensionless quantities that are independent of the system size, apart from 
finite-size corrections. The wrapping probability is an example. It is defined in graphical rep¬ 
resentations as the probability that there exists a connected component which wraps around the 
periodic boundaries, reflecting topological properties of the system under study. Various wrap¬ 
ping probabilities have been studied for the percolation model [1, 2, 3, 4, 5, 6, 7], the Potts 
model [8, 9, 10], loop models [11] etc. They can be related to emergent states of matter such as 
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superfluidity - the superfluid density can be calculated using the winding number [12]. Many di¬ 
mensionless quantities are universal at the critical point, in the sense that their critical values are 
independent of the short-range interactions, lattice types etc. Though their values still depend on 
factors such as the system shape and boundary conditions, the universal dimensionless quantities 
can be used to characterize universality classes of continuous phase transitions [13, 14]. The uni¬ 
versality of the wrapping probability has been proven to depend only on the aspect ratio and the 
boundary twist of the system [1, 2, 3]. Wrapping probabilities are also very useful in determining 
the critical temperature for continuous phase transitions [4, 5]. The well-known Binder ratio [15] 
is another kind of universal dimensionless quantities whose properties have been widely studied 
[16]. 

In a recent work [7], Hu et al. studied the bond (site) percolation model with a fixed total 
number of occupied bonds (sites). The system is referred to as in the canonical ensemble (CE), 
while the model without the constraint is said to be in the grand-canonical ensemble (GCE). 
It is found that, while the wrapping probabilities and dimensionless ratios share the same val¬ 
ues within each of the two different ensembles, some universal quantities in the GCE become 
nonuniversal in the CE, such as the excess cluster number. The percolation model is special, 
since the occupation of each edge is independent of other edges and the average bond density 
has no finite-size dependence. Thus we ask whether the universality of observables depends on 
the ensembles in other models, for which there exist thermal fluctuations that are accompanied 
by finite-size dependence of the particle density. We choose to study the Potts model [17] in the 
random-cluster representation, in which the bonds are treated as particles. Our emphases are on 
the universality and finite-size scaling (ESS) of the wrapping probability in the CE, since their 
universal values at the critical point are exactly known in the GCE [8], and their ESS in the GCE 
is well understood. Moreover, the fact that the wrapping can be defined for each configuration 
significantly simplifies theoretical derivations, as shown later in the text. 

Eor systems in the CE, there exists a constraint that the total number of particles is fixed. 
Erom Eisher renormalization theory [18], under this kind of constraint, the system changes its 
way approaching the critical point as the temperature goes to its critical value, comparing with the 
asymptotic behavior in the GCE. There is a relation between the temperature field in the CE and 
that in the GCE. It turns out that there may exist universal changes for critical exponents in such 
constrained systems: if the specific heat is divergent in the GCE, namely C <x.\T — Tc\^‘‘ with a > 
0, it will become finite in the CE, with = -a/{I - a), where “can” indicates quantities in the 

CE. Other critical exponents may also be renormalized when a > 0, e.g. yScan = y6/(l - a), 7am - 
y/( 1 -a), Vcan - vl{\ —a), where j3, j and v are exponents for the magnetization, susceptibility and 
correlation length, respectively. Eisher renormalization has been extensively studied and applied 
to explain experimental results for constrained systems (see e.g. [7, 14, 19, 20, 21, 22, 23, 24] 
and references therein). However, the study of the universality of dimensionless observables in 
the CE appears to be largely neglected in the literature [19]. 

We conduct a ESS analysis for the Potts model, and support our analysis by extensive numer¬ 
ical simulations. Efficient cluster Monte Carlo algorithms are designed for the simulations in the 
CE . The article is organized as follows. In Sec. 2, we introduce the model and observables. We 
apply Eisher renormalization to finite-size systems in Sec. 3. The universal change of the wrap¬ 
ping probability, and new finite-size corrections are derived in Sec. 4. Sec. 5 presents simulation 
results that verify our derivations. A brief summary and discussion is given in Sec. 6. 


2 


2. Model and Observables 


We introduce the Potts model [17] and define the observables in this section. 

For a lattice G with edge set {e,j), the ^-state Potts model is defined by the Hamiltonian 

9ilkBT ( 1 ) 

eij 

where cr, takes one of the colors 1,2, • • • ,q, representing a spin at site i. Ferromagnetic couplings 
K > 0 occur between nearest-neighbor Potts spins. Performing the Kasteleyn-Fortuin (KF) 
mapping [25], which puts a bond on each edge eq with probability p = (1 - ons gets 

the partition sum of the Potts model in the random-cluster (RC) representation as 

Zrc(v, {v ^ - 1), (2) 

JIQC 


where the summation is over all subgraphs ^ of G. A cluster on a subgraph is defined as a 
connected component consisting of sites and bonds. Afb and JVc are the number of occupied bonds 
and clusters, respectively. The parameter q can take real positive values in the RC representation. 
In two dimensions, the model undergoes a continuous phase transition for q < 4. The critical 
point on the square lattice is predicted to be Vc - -\/q, with the critical bond density pc - 1/2. 
The occupation of different edges are correlated generally for ^ 1. In the limit ^ ^ 1, the Potts 

model reduces to the bond percolation model, for which the presence of bonds at different edges 
are independent. 

In the RC representation, we study the wrapping probabilities R^, Ri, Rh, Re'- Rx counts the 
probability that there exists a cluster connecting to itself along the x direction, irrespective of 
the y direction; Ri is for along the x direction, and with no cluster connecting to itself along 
the y direction; Ri, is for simultaneous wrapping in both directions; and Rg is for wrapping in 
at least one direction. Since Re = 2Ri -i- Rh and - Ri + Rb, only two of these probabilities 
are independent. Exact results for these wrapping probabilities can be obtained from Ref. [8] for 
Potts model in the GCE. 

The relation between wrapping probabilities in the CE and the GCE is presented as follows. 
Denoting R = 1 (0) for a subgraph with (without) a wrapping cluster, we define the average of a 
wrapping probability R^^^ for a finite lattice with linear extension L in the CE as 


L) 




(3) 


where the summation is on all subgraphs which have TVb = pNe bonds. Here, p is the bond 
density, and Ne the total number of edges. The wrapping probability in the GCE is defined as 


R(v, L) 




(4) 


With the definition 


F(p', V, L) 




(5) 
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it can be obtained from Eqs. (3) and (4) that 

R(v, L) = J dp/?“"(p, L)F(p- V, L) , (6) 

where the integral denotes summation over all possible values of p. The quantity F(p; v', L)dp 
is the grand-canonical probability that, for given parameters (v', L), a configuration has a bond 
density p. 

3. Fisher renormalization in finite-size systems 

3.L Finite-size scaling in the grand-canonical ensemble 

For the Potts model in the GCE, it is expected [26] that near the critical point Vc the grand 
potential scales as 

Q(v, L) Q,(v) H- L~‘‘D.s{tU'), (7) 

where Qr and Qj are the regular and singular parts of the potential, respectively. The function Qj 
is universal. The parameter t is the thermal scaling field, and y, = 1/v is the associated exponent. 
One has t ^ A{v — Vc), with A being a nonuniversal metric factor. The magnetic scaling field and 
irrelevant scaling fields are omitted for simplicity. 

The average bond density is 


p(v, L) = ^ + ^lyt ^ (8) 

ov 

where ps{tU') - —Q.\{tU'), pr{v) - -Q'(v'). And its fluctuation scales as 

<(p - pf) - L-\p'Xv) + A^L^y-‘^p[{tLy)) . (9) 

The wrapping probability scales as 

R{v,L)=:^ Ritiy) , (10) 

where ^ is a universal function. 

With Pc = p(vc, L oa) = Pr{vc), near v - Vc, one can define 

Ap = p(v, L)-pc^ (p'c + A^L^y^‘‘p',)Av -I- Aiy-'^ps , (11) 


where Av = v - Vc, p'^ - p[(vc), p's - Ps(0) and ps - Ps(0). This is obtained by Taylor expanding 
Eq. (8) around v'c and ignoring non-linear terms. 

Substituting Eq. (11) into t ^ AAv, one derives that 


AAp - A^U'-^ps 

Pc -t- A^L?-y'~‘^p'^ 


( 12 ) 


In the two-dimensional q - 2 (Ising) case one has 2y, - c/ = 0, and the system in the GCE has 
a logarithmic specific-heat anomaly which scales as Cz,=oo(f) ~ I In |f|| [27]. At the critical point, 
Eq. (9) is replaced by 

<(p-p)2)-L-‘'(Co + CilnL), (13) 
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where Co, Ci are nonuniversal constants. We note that, from Eq. (9) to Eq. (13), Ci contains a 
contribution from p', thus the logarithmic factor InL is associated with p' = -Q"(0), i.e. the 
second derivative of the singular part of the free energy. Accordingly, Eq. (11) changes to 

Ap (Co + Cl In L)Av + AU'~‘‘p, . (14) 


It follows that 

AAp - A^U'^'^ps 
Co + Cl In L 


(15) 


3.2. Finite-size scaling in the canonical ensemble 

In the CE, we conjecture that near the critical density pc, the Helmholtz free energy scales as 

fip, L) - Up) + L-‘^Mx(p, L)). (16) 

In principle, / is related to the grand potential by a Legendre transform / = Q. + vp, with 
p = p(v, L). Thus we suppose that in the CE, x(p, L) is given by tV"'. 

Eor 2y,-d> 0, with p^ + A^U’-'^p^ ^ A^Uy<-‘‘p'^ as L becomes very large, it can be derived 
fromEq. (12) that, 

x(p, L) = tiy - ApL‘‘-ylAp', -pjp', = - 6 , (i?) 

where t = BAp with B - 1/Ap', yr - d - y,, and d = Psip's- The parameter d is universal, 
since both ps and p' are universal quantities. Erom Eq. (17), as L —> oo, r = 0 corresponds to 
t - 0, but tU' - 0 leads to = 6. We note that the hyperscaling relations a - 2 - d jyt and 
Q'can = 2 - djyr, when combined with y^ - d — yt, reproduce the Eisher renormalization result 


ffcan = -a/{I - a). 

Thus, the wrapping probability should scale as 

/;“"(p, L) ^ -S), (18) 

which at p = p(v, L) leads to 

5^f;"“(p, L)ldp^ ^ B^L^y^FF’^'^'\TLy- - (5) = - 6)IA^pf. (19) 

Eor 2y, - d < 0, from Eq. (12), one gets 

x(p, L) = tiy- ^ AApRy- lp[ = , (20) 

where r = BAp with B - Ajp'^, yj - y,. The wrapping probability takes the same form as Eq. 
(18) with d = 0, and Eq. (19) is replaced by 

«9^/?“"(p, L)/5p^ ^ {Tiy-) = . (21) 

Eor 2y, - d - 0, from Eq. (15), one has 

x(p, L) = tiy- AApiy- /Cl In L = / In L , (22) 


where t = BAp with B - A/Ci, yj - yt- Thus the singular part of the free energy changes to 
fsiTL^y ! \TiL), and the wrapping probability scales as 

f;“"(p, L) ^ .R“"(tL"V In L) , 
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( 23 ) 



leading to 


L)/5p^ = L^y'{\nLy^R^'^'^"{ tU^ j \nL)A^ jC] . (24) 

We also note that, close to pc, as L ^ oo, one has |ln|Ap|| oc InL. Thus with t = BAp, 
Eq. (22) leads to 


|fK|T||ln|T|r‘, (25) 

which is the Fisher renormalization result in the thermodynamic limit [18, 20] for the logarithmic 
case. 


4. Universal difference and finite-size corrections 


Expanding the right-hand side of Eq. (6) near p = p(v, L), one gets 

^2Dcan/^ T\ 

Riv, L) - f?^“(p, L) + \ {p - pf), (26) 

where the first order term vanishes since (p -p) = 0. We shall make use of this equation to study 
the difference between critical wrapping probabilities in the CE and the GCE, which are defined 
as s f?‘^™(pc, L ^ oo) and Rq = R{vc, L oo), respectively. 


4.1. 2y,-d>0 

For finite-size systems at the critical point Vc, one has tU' - 0, and from Eq. (17) tL'' ^ 6. 
From Eqs. (10) and (18), one obtains R{vc,L) ^ R{Q) and R‘^’^''(pc, L) ^ .R‘^“(-d), where the 
latter value is different from /?‘^‘“’(p(vc, L), L) ^ .R‘^™(0). In the limit L ^ oa, these approximate 
relations become exact, thus one gets Rq = R(Q) and = .R‘^™(-d). 

The combination of Eqs. (9) and (19) at the critical point Vc leads to 

52p“n(p,L) _ .R“""(0) L‘^^2y,^can"(0)p, 

2dpi > 2p; ^ 2A2pf ■ 

Substituting the above results into Eq. (26), we derive 


Rq ^ .R“"(0) -H 


^can''(0) 

2p' 


L'''2>Vpcan"(Q)p/ 

2A2pf 


(28) 


Equation (28) tells that: (i) the critical values of the wrapping probability in the GCE are 
different from the values in the CE. This remains true in the limit of L ^ oo, i.e. Rq + = 

pcan(_^) since Rq, p' and 5 are universal, (ii) .R‘^™(0) and .^“^“"(O) are universal, as well as 
the difference between 7?^“" and Rq. It can be derived that higher-order terms omitted in Eq. (26) 
lead to universal terms similar to (0)/2p', but with high-order derivatives of These 
terms also contribute to the difference between Rq and 

The leading finite-size correction in the GCE has an exponent y,, i.e. a term proportional to 
V' is present on the left-hand side of Eq. (28). Since one usually has yi + d - 2y,, there should 
exist a correction term ~ L''™ in the CE with jcan - d - 2y,, which cancels the term proportional 
to in Eq. (28). 
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4.2. 2yt — d < 0 

For 2y, — d < 0, at the critical point Vc, one has tL'' = 0. Following a procedure similar to 
that for 2yt - d > 0, we get 


Ro - Rr + 




r 2y,-d 


(29) 


Thus, at criticality we have Rq - the wrapping probability is identical in the two ensembles, 
and a finite-size correction with exponent jcan = 2y, - c/ occurs in the CE. 


4.3. 2y,-d^0 

For 2y, - d - 0, one has j In L = 0 at the critical point Vc. It can be derived that 


Rq 


Dcan 

Rq ■ 


^can 

2Ci 


/InL 


(30) 


which tells that Rq - and a logarithmic correction term proportional to 1 / In L is present in 
the CE. 


4.4. Remark 

The universal difference between the wrapping probability in the CE and GCE for systems 
with 2yt - d > 0 is the main finding of this work. It is important to note that, for this universal 
difference, there is a contribution coming from fluctuations of the bond density in the GCE, even 
when the universal parameter S is zero. The suppression of bond-density fluctuations also leads 
to finite-size corrections with an exponent jcan = -|2yf - d\ for systems with 2y, — d + 0, and 
logarithmic corrections such as a term ~ 1 / In L for 2y, — d - Q. Einite-size corrections can 
also come from Eisher renormalization [22, 23, 24]: for 2yt - d > 0, from Eq. (12) to (17), a 
correction term with should be present when higher-order effects are considered; similarly 
one has a correction term ~ for 2y, — d < 0, and ~ 1 / In L for 2y, — d - 0. 

We mention that the above findings should apply to other universal observables like Binder 
ratio and to other canonical-ensemble systems like the dilute Potts model with a fixed total num¬ 
ber of vacancies. 


5. Numerical verification 

To verify the analysis in the previous section, we conducted extensive Monte Carlo calcula¬ 
tions, for which the simulation details and results are presented below. 

5.1. Monte Carlo Algorithms 

Eor the simulation of the Potts model in the GCE, we employ the Swendsen-Wang (SW) 
algorithm [28]. In the CE, the total number of bonds is fixed and we designed the following 
algorithm for updating the configurations: (i) randomly distribute the bonds on the lattice edges; 

(ii) for each cluster, color it to be ‘k’ with probability pk - Ijq for each color {k - 1,2,..., q)-, 

(iii) independently on each subgraph G{Vk\ (Vk represents all sites with color ‘k’, and G[Vk] 
consists of Vk and the edges between sites with color ‘k’), employ Kawasaki dynamics [29] for 
bond percolation, i.e. exchange the state of two randomly selected edges; (iv) erase the colors; 
(v) repeat steps (ii), (iii) and (iv) until the required number of samples is reached. The canonical 
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Table 1: Details of the simulations for the Potts models at the critical point. The number of updates between two samples 
were chosen so that the samples are approximately independent. For each size, the total number of samples were taken 
from over 50 independent jobs, for which each job had lots of updates (roughly equivalent to 1/5 of the number of 
samples) thrown for thermalization before sampling. 


Model 

L 

Number of samples 

q^3 GCE 

(8,16,32,64,128,256,512) 

3.0 X 10** 

^ = 3CE 

[8,16,32,64,128] 

1.2 X 10® 


[256] 

2.6 X 10** 


[512] 

1.4 X 10** 


[1024] 

4.0 X 10® 


[2048] 

2.3 X 10® 

q^2 GCE 

[8,16,32,64,128,256] 

5.0 X 10® 


{512] 

2.5 X 10® 

^ = 2CE 

[8,16,32,64,128,256] 

5.0 X 10® 


[512] 

2.5 X 10® 


algorithm can be easily adapted to simulate RC models with non-integer q value {q> 1): for each 
cluster, color it to be ‘1’ with probability p - Ijq, and ‘2’ with p - \ — Ijq-, then on subgraph 
G[Vi], employ Kawasaki dynamics for bond percolation. For non-integer q'-state RC model in 
the GCE, the Swendsen-Wang-Chayes-Machta algorithm (SWCM) [30] algorithm can be used. 

5.2. Simulation details 

The simulations were conducted on Lx L square lattices with periodic boundary conditions. 
We simulated the q = 3 Potts model which has ly, - d = 2/5 > 0, and the q - 2 Potts model 
which has 2y, — d - 0. For the models in the GCE, the critical point is v'c = yfq, and the critical 
bond density is p* = (7Vb)/2L^ = 1/2. Eurther details of the simulations at the critical point are 
summarized in Table 1. With L - 64,128,256, we also did simulations for the models at several 
points near the critical point, for which over 10^ samples were taken at each point for each size. 

5.3. Numerical results for 2y, — d > 0 

In the CE, the thermal renormalization exponent is renormalized as = d - yt - 4/5 for 
q - 3 Potts model. Prom the simulation data, we plotted this scaling renormalization for two 
independent wrapping probabilities and Rh, as shown in Pig. 1. The figure also shows that the 
critical wrapping probabilities in the CE are different from those in the GCE. 

Using the method of least squares, we performed fits with our Monte Carlo data at the critical 
point by the formula 

R^Ro + b^U' +b2U\ (31) 

where Rq is the universal value of the wrapping probability; yi andy 2 are the leading and sublead¬ 
ing correction exponents, both being negative; and bi, b 2 are nonuniversal amplitudes. Por the 
GCE, the leading correction exponent is yi - y, = -4/5 [31, 32]. And for the CE, the correction 
exponents are expected to be yi = yean = -2/5 and y 2 = -4/5. We set a lower cutoff L > 
on the data and observed the change of the residual as Lmin increases. Subsequent increases 
of Lmin do not make the value drop vastly by more than one unit per degree of freedom. 
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Figure 1: Wrapping probabilities Ri, (top), (bottom) versus L‘^~^'(p — Pc) in the canonical ensemble (left), and versus 
U‘(v - Vc) in the grand-canonical ensemble (right), for the q = "i Potts model on square lattices, with y, = 6/5 and the 
spatial dimension d = 2. Data points for three different sizes L = 64,128,256 (cross, up-triangle, down-triangle) collapse 
well into a single curve. Solid lines indicate the coordinates of the critical values. 


Table 2 summarizes the results of our fits. The fits were made by assuming one or two 
correction terms, and in alternate lines using the predicted values of y\ and y 2 - For the GCE a 
fit was made with Rq fixed at its theoretical value. As expected from the analysis in Sec. 4 , the 
critical values of the wrapping probabilities in the CE are significantly different from those in the 
GCE. It is interesting to note that the predicted leading correction amplitude bi is consistent with 
bi - 0, while this no longer holds for the dilute q=3 Potts model (see the manuscript later). 


5.4. Numerical results for 2y, — d — 0 

Eor the logarithmic case 2y,-d - 0, near the critical point, the analysis in Sec. 3 tells that the 
wrapping probabilities scale as U'(p - pc)/ InL in the CE, in contrast to U'{v - Vc) in the GCE. 
This scaling is shown in Eig. 2 for the Ising {q -2 Potts) model on the square lattice. 

Eor the Ising model at the critical point, we performed fits to our data in the GCE with the 
formula 

R^Rq + biU' , (32) 

and in the CE with the formula 


„ ^1 ^2 bi 

R — Ro -t---F- — + - — . 

InL (lnL)2 (InL)^ 


(33) 


The fit results are summarized in Table 3. The value of Rq was obtained from the fit or fixed at 
its theoretical value in alternative lines, and for the CE only two correction terms were assumed 
when Rq was not fixed. It can be seen that the wrapping probabilities share the same universal 
values in the CE and the GCE, and that the leading finite-size correction term is proportional to 
1 / In L in the CE. Eigure 3 demonstrates the logarithmic correction by showing the data for Ri. 
These results are consistent with the analysis Sec. 4. 
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Table 2: Fit results for wrapping probabilities of the ^ = 3 Potts model. The abbreviations “GCE” and “CE” stand for the 
grand-canonical and canonical ensemble, respectively. Entries indicate that the con'esponding term is not included 
in the fit, and the numbers without error bars are fixed in the fits. Exact values for Rq were obtained using formulas in 
Ref. [8]. EiTor margins are quoted as two times of the statistical errors in the fits. 





Eits using R - 

Ro + bxU' 

+ b2U^- 






Ro 

b\ 

y\ 

bi 

y2 

^min 

x^Idf 

Rx 

GCE 

0.695 000176 

-0.04(2) 

-1.07(18) 

- 

- 

16 

0.2 



0.695 2(3) 

-0.021(12) 

-4/5 

- 

- 

64 

0.1 


CE 

0.667 2(3) 

-0.05(2) 

-0.69(13) 

- 

- 

64 

0.2 



0.667 3(5) 

-0.005 (6) 

-2/5 

-0.05(2) 

-4/5 

64 

0.1 

Ri 

GCE 

0.118 666 330 

0.056 (6) 

-0.86(3) 

- 

- 

32 

0.6 



0.11868(10) 

0.040(6) 

-4/5 

- 

- 

128 

0.5 


CE 

0.18016(14) 

0.152(14) 

-0.77(3) 

- 

- 

64 

0.8 



0.1801(2) 

0.003 (3) 

-2/5 

0.157(10) 

-4/5 

64 

0.7 

Rh 

GCE 

0.576 333 845 

-0.115(18) 

-0.97(6) 

- 

- 

16 

0.4 



0.576 6(2) 

-0.066 (6) 

-4/5 

- 

- 

32 

0.1 


CE 

0.487 0(4) 

-0.20(4) 

-0.75(5) 

- 

- 

64 

0.3 



0.487 2(6) 

-0.008 (8) 

-2/5 

-0.21(3) 

-4/5 

64 

0.2 

Re 

GCE 

0.813 666 506 

0.034(9) 

-0.84(9) 

- 

- 

16 

0.7 



0.813 8(2) 

0.027 (4) 

-4/5 

- 

- 

32 

0.3 


CE 

0.847 30(14) 

0.103(10) 

-0.82(4) 

- 

- 

32 

0.2 



0.8474(3) 

-0.001 3(28) 

-2/5 

0.102(7) 

-4/5 

32 

0.2 



0.54 

0.5 

0.46 

0.42 

0.68 

0.64 

0.6 

0.56 


- 0.4 - 0.2 0 0.2 - 0.2 0 0.2 0.4 

L(p-p^)/lnL L(v-v^) 


Figure 2: Wrapping probabilities Rt, (top), Rx (bottom) versus - pc)I InL in the canonical ensemble (left), and 

versus L'^(v - Vc) in the grand-canonical ensemble (right), for the Ising model on the square lattice, with y, = 1 and 
the spatial dimension d = 2. Data points for three dilferent sizes L = 64,128,256 (cross, up-triangle, down-triangle) 
collapse well into a single curve. Solid lines indicate the coordinates of the critical values. At p = pc, a small deviation 
of from its universal value 0.480 comes from finite-size corrections, for which the leading term scales as 1/ InL. 
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Table 3: Fit results for the wrapping probabilities of the Ising model. The abbreviations “GCE” and “CE” stand for the 
grand-canonical and canonical ensemble, respectively. Entries indicate that the con'esponding term is not included 
in the fit, and the numbers without error bars are fixed in the fits. Exact values for Rq were obtained using formulas in 
Ref. [8]. Error margins are quoted as two times of the statistical errors in the fits. 




GCE, Fits using R - Rq + b\U^ 




Ro 

b\ 



^min 

X^IDF 

Rx 

0.627 18(15) 

-0.13(21) 

-2.1(8) 


8 

1.1 


0.627138 794 

-0.18(26) 

-2.2(7) 


8 

0.9 

Ri 

0.14643(5) 

0.166(20) 

-1.60(6) 


8 

0.6 


0.146436 927 

0.168(16) 

-1.60(4) 


8 

0.5 

Rb 

0.48075(18) 

-0.26(9) 

-1.68(16) 


8 

1.1 


0.480701867 

-0.27(8) 

-1.71(12) 


8 

0.9 

Re 

0.773 62(17) 

0.09(4) 

-1.47(24) 


8 

1.3 


0.773 575 721 

0.08(3) 

-1.43(16) 


8 

1.1 


CE, Fits using R - Rq + bi/ \nL + b 2 /(lnLY + b^j (In LY 



Ro 

b\ 

bi 

^3 

^min 

xVdf 

Rx 

0.626 1(12) 

-0.030(9) 

-0.013(16) 

- 

16 

0.4 


0.627138 794 

-0.042(3) 

0.029(14) 

-0.05(2) 

8 

0.3 

Ri 

0.148 6(24) 

0.134(24) 

0.02(6) 

- 

64 

2.1 


0.146436 927 

0.162 1(24) 

-0.104(18) 

0.18(3) 

16 

1.0 

Rh 

0.477 9(28) 

-0.166(25) 

-0.03(6) 

- 

32 

0.8 


0.480701867 

-0.205(6) 

0.14(5) 

-0.25(8) 

16 

0.5 

Re 

0.775 6(23) 

0.095(21) 

0.03(5) 

- 

32 

0.1 


0.773 575 721 

0.118(3) 

-0.053(14) 

0.097(20) 

8 

0.3 



Figure 3: Quantity l/(/?i - Ri,o) versus InL for the Ising model in the canonical ensemble. The solid straight line is 
described by !/(/?] - ^i,o) = 6.70 InL, with Ri^q = 0.146 436 927 being the exact value obtained using formulas in 
Ref.[8]. Small deviation of data points from the straight line can be attributed to higher-order correction terms. 
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5.5. 2y, — d <0 

For 2yt - d < 0, an example in the CE is the percolation model with fixed number of bonds 
(sites). Results for two-dimensional percolation models in Ref. [7] are consistent with our an¬ 
alytical results for systems 2y, - d < 0, i.e. values of the universal wrapping probabilities and 
dimensionless ratios in the CE are the same with those in the GCE; and correction exponent 
Jean = -|2y/ - d\ is present in the ESS. 

5.6. The dilute Potts model 

We also investigated the q - 3 dilute Potts model to demonstrate the universality of the 
wrapping probabilities in the CE. The dilute Potts model is defined by the Hamiltonian 


'HdlkeT = 2 - 5^,.o) - ^ . 


(34) 


Here cr, = 0 stands for a vacancy at site i, and cr, = 1,2,...,^ for one of the Potts states. The 
abundance of vacancies is controlled by the chemical potential D. In the limit D -oo, the 
vacancies are excluded and the above Hamiltonian describes the pure Potts model. The partition 
sum of the dilute Potts model in the RC representation is 



(35) 


In comparison with the RC representation of the pure Potts model, a term appears, with 
jj. - e^, and a number TVy of vacancies. One has TVy = 0 for the pure Potts model. In two 
dimensions, for q < A, the phase diagram of the dilute Potts model in the (v, fi) plane consists 
of a line of Potts critical transitions and a line of first-order transitions, which join at a tricritical 
point [33]. As a single fixed point governs the Potts critical transitions, the dilute Potts model 
belongs to the same universality class as the pure Potts model. We define the dilute Potts model 
with a fixed number of Potts spins (the number of vacancies is also fixed since the total number 
of sites is conserved) as the model in the CE. Our finite-size analysis for the pure Potts model 
applies also to the dilute Potts model, with the bond density replaced by the vacancy density, and 
the thermal scaling field being t ^ AQi-fic) {t should be approximated by Ai(ju - fic) + A 2 {v - Vc), 
but here we take v = )■ 

To simulate the dilute Potts model, in the GCE, the SW method is used for updating Potts 
spins, and the Metropolis algorithm is employed to allow fluctuations between the vacancies and 
Potts spins. In the CE, the S W method is also used for updating Potts spins, but in order to fix the 
total number of vacancies while still allow for their spatial fluctuations, a Kawasaki-like algo¬ 
rithm is conducted, which exchange the states of two randomly selected sites with a probability 
satisfying the detailed balance condition. Replacing the SW method by the SWCM method, this 
algorithm for the CE applies generally to dilute RC models with q > 1. Eor integer q'-state di¬ 
lute Potts models in the CE, the geometric cluster algorithm [34] can be used, which has a more 
efficient dynamics than the Kawasaki-like algorithm described above. The simulation results 
at the critical vacancy density in the CE were obtained from linear interpolations between 
Mm = [Pc^^] and Afy = [Pc^^] + where [■] denotes the integer part of the number. 

Numerical simulations were done at a critical point Kc = 1.169 41(2), Dc - 1.376 483 (5) 
[35] for the q - 3 dilute Potts model, with the average vacancy density being p^ = 0.105 28(1) 
[35]. This point was determined such as to suppress the leading irrelevant scaling field. Thus, in 
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Table 4: Fit results for wrapping probabilities of the q = 3 dilute Potts model. The abbreviations “GCE” and “CE” 
stand for the grand-canonical and canonical ensemble, respectively. Entries indicate that the corresponding term is 
not included in the ht, and the numbers without eiTor bars are fixed in the fits. Exact values for Rq were obtained using 
formulas in Ref. [8]. Error margins are quoted as two times of the statistical eiTors in the fits. 


Fits using R - Rq + b\U' + b 2 U- 




Ro 

b\ 

yi 

b2 

y2 

^min 

X^DF 

Rx 

GCE 

0.695 000176 

-0.16(8) 

-2.1(3) 

- 

- 

8 

0.6 



0.695 04(6) 

-0.121(6) 

-2 

- 

- 

8 

0.6 


CE 

0.667 7(12) 

0.045(15) 

-2/5 

-0.26(4) 

-4/5 

32 

0.5 

Ri 

GCE 

0.118 666 330 

0.29(12) 

-2.11(15) 

- 

- 

16 

0.9 



0.118 67(3) 

0.19(4) 

-2 

- 

- 

16 

1.1 


CE 

0.179 7(3) 

-0.101(3) 

-2/5 

0.153(7) 

-4/5 

16 

0.5 

Rh 

GCE 

0.576 333 845 

-0.40(8) 

-2.08(9) 

- 

- 

8 

0.6 



0.576 39(7) 

-0.342(6) 

-2 

- 

- 

8 

0.7 


CE 

0.487 9(14) 

0.15(2) 

-2/5 

-0.42(5) 

-4/5 

32 

0.5 

Re 

GCE 

0.813 666 506 

0.09(4) 

-1.95(23) 

- 

- 

8 

0.6 



0.813 68(5) 

0.099 (5) 

-2 

- 

- 

8 

0.6 


CE 

0.847 6(10) 

-0.058(13) 

-2/5 

-0.1(4) 

-4/5 

32 

0.4 


the GCE, we expect corrections with a scaling exponent smaller than y, = -4/5, e.g. the integer 
-2 [35]. We conducted the simulations at 7 different sizes with 8 < L < 512. The number of 
samples were about 1.7 x 10^ for L < 64, 2.0 x 10* for L - 128,256, and 6.4 x 10^ for L = 512 
in the GCE; and about 5x10^ for each size in the CE. 

We performed fits for the data of the q - h dilute Potts model using ansatzes that are similar 
to those for the q — 2> pure Potts model. The fit results for the wrapping probabilities are shown in 
Table 4. Similar to the case for the pure q - 3 Potts model, the data in the GCE can also be well 
described by known exact values of the wrapping probabilities, as expected from universality; 
and the wrapping values in the CE are different from those in the GCE for the q - 3 dilute Potts 
model. Moreover, in the CE, within uncertainties, the wrapping values for the q = 3 pure and 
dilute Potts model agree with each other. This demonstrates that universality of the wrapping 
probability still holds in the CE. In contrast to the small finite-size corrections in the GCE, the 
data in the CE are consistent with the appearance of correction terms with exponents -2/5 and 
-4/5, as expected from the analysis in Sec. 4. As an example, Eig. 4 illustrates the finite-size 
behavior of Ri. 


5.7. Dimensionless ratios 

We also observed two dimensionless ratios Qc, Qs which are defined as following; 

(Cl >2 


Qc^ 


{Ch 


(36) 


where Ci is the size of the largest cluster; 


Qs 


{q-\){S2? 

{{q+DSl-lS^y 

13 


(37) 









Figure 4: Wrapping probability R\ versus for the dilute q = 3 Potts model in the canonical ensemble (main plot) 
and the grand-canonical ensemble (inset), where it is assumed that the leading correction exponent y\ takes value -2/5 
and -2, respectively. The (blue in the web version) triangular point represents the critical value of Ri in the canonical 
ensemble = 0.179 7(3). The solid line connecting data points in the canonical ensemble is described by R\{L) = 
0.179 7 - O.lOlL-2/5 + 0.!53L^''/^ and the line in the inset is Ri(L) = 0.118 666 330 + 0.21L-^. 


where Sk - Y/jh with sj being the size of the jth cluster divided by the volume For 
q -2, the latter quantity is equivalent to the Binder ratio [15]. It is found that these ratios behave 
similar to the wrapping probability. For both the q - pure and dilute Potts model, which have 
2y, - d > 0, the critical values of the ratios in the CE are different from those in the GCE, but 
they remain universal. Table 5 summarizes the critical values of the two dimensionless ratios and 
the wrapping probabilities for the q - 3 Potts models. 

6. Summary and discussion 

In this work, based on ESS analysis, we derive that for systems with 2y, - > 0 in the GCE, 

in the limit of L —> oo, the universal wrapping probabilities at criticality become different in the 
CE. However, they are still universal in the CE. Eor 2yt - d < 0, the critical universal wrapping 
probabilities do not change. It is also derived that, for 2y, — d + 0, finite-size corrections with 
exponent Jean = —\2yt—d\ appear in the CE. Eor 2yt-d = 0, i.e. the specific heat has a logarithmic 
anomaly in the GCE, the leading correction term changes to a logarithmic form in the CE, such as 
1/ InL for the two-dimensional Ising model. Other dimensionless quantities, such as the Binder 
ratio, behave similar to the wrapping probabilities. These results are supported by numerical 
results from extensive Monte Carlo simulations. 

The simulations were conducted for two-dimensional critical systems, including the q — 3 
Potts model {y, — 6/5), the Ising model {y, - 1), and the percolation model iy, — \jv - 3/4). 
As described in Sec. 2, these models can be described in the RC representation of the ^-state 
Potts model, which are well defined also for positive noninteger q. In two dimensions, exact 
critical exponents for the Potts model can be obtained from the Coulomb gas theory [31] or the 
conformal invariance theory [32]. In the former case, they can be expressed in the Coulomb gas 
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Table 5: Critical value of the wrapping probabilities and dimensionless ratios in the thermodynamic limit L oo, for 
the q = 3 pure and dilute Potts model in the grand canonical ensemble (GCE) and the canonical ensemble (CE). These 
results tell that the CE values are different from the GCE values, and the differences ai‘e universal. The exact values were 
obtained using fonnulas in Ref. [8]. 


model 

Rx 

Ri 

Rb 

GCE pure 

0.695 2(3) 

0.118 68(10) 

0.576 6(2) 

GCE dilute 

0.695 04(6) 

0.118 67(3) 

0.576 39(7) 

GCE exact 

0.695 000 176 

0.118 666 330 

0.576 333 845 

CE pure 

0.667 3(5) 

0.180 1(2) 

0.487 1(7) 

CE dilute 

0.667 7(12) 

0.179 7(3) 

0.487 9(14) 

model 

Re 

Qc 

Qs 

GCE pure 

0.813 8(2) 

0.936 36(14) 

0.854 3(3) 

GCE dilute 

0.813 68(5) 

0.936 26(14) 

0.854 1(3) 

GCE exact 

0.813 666 506 



CE pure 

0.847 4(3) 

0.974 2(2) 

0.919 2(5) 

CE dilute 

0.847 6(10) 

0.974 3(4) 

0.919 5(9) 


coupling constant g which depends on q by 

^ = 2 + 2cosy, (38) 

with 2 < g < 4 for the critical Potts model. The thermal renormalization exponent y, is given by 

X,^d-y,^--l . (39) 

8 

Thus, for the critical Potts model in two dimensions, one has 2y,-d > 0 for 2 < ^ < 4, 2y,-c/ = 0 
for^ = 2, and2y,-^f < OforO < q <2. In higher dimensions, the Ising model has2y^-^f ^ 0.174 
for d - 3 [36], and the percolation model has 2y, - d <0 also for > 2 [5, 37]. The results 
summarized in the previous paragraph should apply but not restricted to these models, as long 
as the wrapping probability can be properly defined (wrapping is trivial on a complete graph). 
It should be mentioned that, for the percolation model, since the particle density has no finite- 
size dependence, the derivation presented need some modifications [7], which do not change the 
results for the wrapping probability. 

Exact values for the critical universal wrapping probability of the Potts model in the GCE 
were obtained through the analysis of the homology group of the torus based on a method in¬ 
troduced by di Erancesco et al. [38, 2, 8, 9, 10] (Ref. [10] contains a brief review). It remains 
an open question whether exact results can be obtained for systems with 2y, - > 0 in the CE. 

Although critical universal wrapping probabilities and dimensionless ratios remain universal in 
the CE, some universal parameters may become nonuniversal, such as the excess cluster number 
for percolation [7], which is universal in the GCE. It would be interesting to study properties of 
other quantities, and constraints other than the constraint on the number of particles, e.g. fixing 
the magnetization. 

Einally, we note that the wrapping probabilities and dimensionless ratios considered in this 
work are different from the universal quantities studied by Izmailian and Kenna [19], which are 
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expressed as combinations of universal amplitudes governing the behavior of various quantities 
near the critical point. 
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